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ABSTRACT 

Dark Stars (DS) may constitute the first phase of stellar evolution, powered by dark 
matter (DM) annihilation. We will investigate here the properties of DS assuming 
the DM particle has the required properties to explain the excess positron and elec- 
tron signals in the cosmic rays detected by the PAMELA and FERMI satellites. Any 
possible DM interpretation of these signals requires exotic DM candidates, with an- 
nihilation cross sections a few orders of magnitude higher than the canonical value 
required for correct thermal relic abundance for Weakly Interacting Dark Matter can- 
didates; additionally in most models the annihilation must be preferentially to lep- 
tons. Secondly, we study the dependence of DS properties on the concentration pa- 
rameter of the initial DM density profile of the halos where the first stars are formed. 
We restrict our study to the DM in the star due to simple (vs. extended) adiabatic 
contraction and minimal (vs. extended) capture; this simple study is sufficient to 
illustrate dependence on the cross section and concentration parameter. Our basic 
results are that the final stellar properties, once the star enters the main sequence, 
are always roughly the same, regardless of the value of boosted annihilation or con- 
centration parameter in the range between c=2 and c=5: stellar mass ~ IOOOM0, 
luminosity ~ 1O 7 L0, lifetime ~ 10 6 yrs (for the minimal DM models considered 
here; additional DM would lead to more massive dark stars). However, the lifetime, 
final mass, and final luminosity of the DS show some dependence on boost factor 
and concentration parameter as discussed in the paper. 
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1 INTRODUCTION 



Spolyar et al.1 (|2008) first considered the effect of DM particles on the first stars during their for- 
mation. The first stars formed when the universe was about 200 million years old, at z = 10 — 50, 
in 1O 6 M haloes consisting of 85% DM and 15% baryons predominantly in the form of H and 
He from big bang nucleosynthesis. The canonical example of particle DM is Weakly Interacting 
Massive Particles (WIMPs). In many theories WEVIPs are their own antiparticles and annihilate 
with themselves wherever the DM density is high. The first stars are particularly good sites for 
annihilation because the y form at high redsh ifts (density scales as (1 + z) 3 ) and in the high density 



centers of DM haloes. 



Spolyar et al.l (|2008|) found that DM annihilation provides a powerful heat 



source in the first stars and suggested that the very first stellar objects might be Dark Stars (DS), 
a new phase of stellar evolution in which the DM - while only a negligible fraction of the star's 
mass - provides the key power source for the star through DM heating. Note that the term 'Dark' 
refers to the power source, not the appearance of the star. Dark Stars are stars made primarily of 
hydrogen and helium with a smattering of dark matter (< 1 % of the mass consists of DM); yet they 
shine due to DM heating. 



In subsequent work, 



Freese et al. 



Spolyar et al. 



(2009) studied the stellar struc- 



(I2008af) and 

ture of the dark stars, and found that these objects grow to be large, puffy (~ 10 A.U.), bright 
(~ 1O 7 L ), massive (they grow to at least ~ 5OOM ), and yet cool (~ 6000 - 10, 000K surface 



temp erature) objects. They last as long as they are fed by dark matter. In another paper, 



Freese et al. 



(|2010h considered the possibility of an extended period of dark matter heating and the consequent 
growth to Supermassive Dark Stars > 1O 5 M . By contrast, in the standard case, when DM heating 

le standard fusion pow ered first stars) U form by accretion 



is not included, Population III stars (t ; 

onto a smaller protostar ~ 1O~ 3 M (lOmukai & Nishilll998|) up to ~ 14OM and surface tem- 



peratures exceeding 5 x 10 5 K. This higher surface temperature in the standard picture inhibits 
the ac cretion of the gas, as various feedback mechanisms become effective at those high tempera- 



tures (McKee & Tan 



2008). For reviews of first stars formation in the standard scenario, where DM 



1 Pop. Ill stars refer to the first stars to form in the universe and are uncontaminated by previous stellar evolution. They consist only of the elements 
produced during Big Bang Nucleosynthesis. 
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heatin g is not included see e.g. 



(2004); 



Ripamonti & Abe l (2005); 



Barkana & Loeb 



Yoshida et al. 



(2001); 



(2006) 



Yoshida et al. 



(2003); 



Bromm & Larson 



WIMP annihilation produces energy at a rate per unit volume 

Qdm = n 2 x (av)m x = (av)p x /m x , (1) 

where n x is the WIMP number density, m x is the WIMP mass, and p x is the WIMP energy den- 
sity. The final annihilation products typically are electrons, photons, and neutrinos. The neutrinos 
escape the star, while the other annihilation products are trapped in the dark star, thermalize with 
the star, and heat it up. The luminosity from the DM heating is 

Ldm ~ f Q f QoudV (2) 



where fq is the fraction of the annihilation energy deposited in the star (not lost to neutrinos) and 
dV is the volume element. 

The DM heating rate in Dark Stars scales as the square of the WIMP density times the anni- 
hilation cross section, as can be seen from Eq.©. The WIMP density inside the star adjusts in 
response to changes in the star's baryonic mass profile, since the gravitational potential well of 
the star is determined by the baryons. (The DM profile responds to changes in the gravitational 
potential due to the conservation of adiabatic invariants.) 

In this paper, we investigate the dependence of Dark Star properties on these two quantities: 
(i) the annihilation cross section and (ii) the density of the halo within which the star forms, as 
characterized by the concentration parameter. 

For a short list of papers by various other authors that have continued the work of ISpolvar et al 



(2008) an d explored t he rep e rcussions of DM he a ting in the first stars see 



(2008); 



Taoso et al 



Ripamonti et al 



cussed in 



(2008) 



(2010); 



Freese et al 



Yoon et all (12 008): 



Sivertsson & Gondolo (2010). 



OOlOt); 



Zackrisson et al 



Ripa monti et al 



(2009) 



Iocco(2008); 



Iocco et al 



Gondolo et al 



toot); 



(2010a 



heir potential observability has been dis- 



bj). The possibility that DM annihilation 



mig ht have effects on today's stars was actual l y considered in the ' 80s by various authors such 



as 



Krauss et al. 



09851); 



Bouquet & Salati 



(198 



Salati &Silkl(ll989h . More recentl y the effect on 



today ' s stars has been re-examined under the assumption that DM is m ade of WIMPs (Moskalenko & Wai 



2007; 



Scott et al. 



2007; 



Bertone & Fairbairn 



inelastic dark matter (IHooper et al 



2010) 



2008 



Scott et al. 



2008) or within the hypothesis of 



© 0000 RAS, MNRAS 000, 000-000 



1.1 Boosted Leptophilic Annihilation Motivated by PAMELA data 



Recent measurements of cosmic ray positrons and electrons in the GeV-TeV range could have 
signi ficant implications on our understanding of dark matter (DM). The PAMELA collabora- 



tion (Adriani et al. 



2009b, 



2010 ) reported a e + flux excess in the cosmic energy spectrum from 
10 to 100 G eV, reinforcing what was previously observed up to an energy of 50 GeV b y the HEAT 



experiment (|Barwick et al 



2009a : b) has found 



19971) . The FERMI-LAT collaboration (lAbdo et all 
an excess in the e + + e~ flux in the 100 — 1000 GeV energy range above the background given 
by a conventional diffusive model, albeit in conflict with a much larger excesses in flux in the 
500 GeV range previously reported by ATIC collaboration @ It is worth mentioning that a simple 
power-law fit of the FERMI-LAT electron energ y spectrum is possible, being consistent w ith both 



astrophysical sources, or DM annihilation (e.g 



Di Bernardo et al. 



2009; 



Grasso et al. 



2009). 



If confirmed, there are several possible explanations for the positron excess. T he signals could 



Profumo 2008 



2009 



Hooper et al 



Mertsch & Sarkar 



2009a; 



200 



9, 



Ahlers et al. 



Shaviv et al. 



2009 



20091) . Uncertainties in cosmic ray propagation in 



or supernova shocks (e.g. 


Boulares 


198* 


Blasi 


2009 


Fujita et al. 


20o4 


Malvshev et al. 



the galaxy leave open the possibility that standard cosmic ray physics could explain the signal 



(|Delahaye et al. 



2009). Another possib 



from DM annihilation (Baltz et al 



2008 



Hooper & Silk 2005 



2009 



2009 



2009 



2009 



Cirelli et al 



Nomura & Thaler 



Hooper & Zure kl2009l:lMeade et alJ2 009 



Bai et al. 



200' 



Chen et al 



Shirai et al. 



200 1 



Bi et al. 



2002 



2008: 



e interpretation of the data cou 



Kane et al 



2002 



Cirelli & Strumi 



Arkani-Hamed et a 



2009; 



Cholis et al 



2009a.b 



2009i 



Ibarra & Tran 



Yin et al. 



2009 



straints against such an interpretation (|Abdo et al 



2009; 



de Boer et al 



2009 



2008 



2009a 



Zurek 2009; 



d be in terms of a signal 



2002; 



Graiek et al 



large r et al 



bk; 



Ibarra et 



Baic et a 



2010; 



201 



0; 



al 



Hooper et al. 



2009J: 



2009 



Cirelli et al 



Meade et al. 



2009; 



2004 



Nelson & Soitzer 



Bergstr om et al. 



Ibe et al. 



Chen et al 



20 



Abazajian et al 



0, 



2009a 



2009; lFox& Poppitz 

T 

2010) or decay (Chen & Takahashi 



Ishiwata et al. 



Others point out con- 



2010). 



The prospect that DM has been detected, although indirectly, has stirred a lot of excitement 
and a flurry of interest. There are several model-building challenges, however. If dark matter is 



2 The ATIC balloon experiment Chang et alj|2008 ) found an increase in the e+ + e 
which has apparently been ruled out by FERMI. 



spectrum between 300 and 800 GeV above background 
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to explain the positron flux excess in the cosmic -ray spectrum, in most models the products of 
decay or annihilation must be primarily leptonic since an excess in the anti-proton fluxes is not 



found (Adriani etal. 



2009ah 



Moreover, in order to explain the observed signals the annihilation rate needs to be on the or- 



ments to the cross-section, such as the Sommerfek 


enhancement ( 


Sommerfeh 


i 193J 


1 

Hisano et al. 


2004 


: Profumo 2005: Marc 


i-Russell et al 


2008|; 


Cirelli et al. 


200S 


); 


Lattan/ 


d & Si] 


llJ 2009 


) or a 


Breit-Wigner enhancement ( 


Feldman et al. 


20091 


Guo & Wu 


20091: 


be et al. 


2009b; 


Kadota et al. 



2010J). S cenarios involving non-standard cosmologies 



tion (e.g 



Gelmini & Gondololl2 008; 



Catena et al. 



lave also bee n considered as a possible solu- 



2009; 



PaUisl2010|) . Alternatively, the annihilation 



rate could be greater than the sta ndard prediction if the dark matter is not smoothly distributed in 



the l ocal halo of the Milky Way (|Diemand et al. 



2008 



Hooper et al. 



2009b; 



Kamionkowski et al. 



2010J). For a review of the DM explanation for the the cosmic ray e ± excess see 



He (2009) 



In light of these new dark matter models, this paper will study the possible mo difications to 



the e volution of a Dark Star (DS) for the case of leptophilic boosted dark matter. In lSpolyar et al. 



(|2009Q a comprehensive study of DS for various WIMP masses was done, but the annihilation 
cross section was kept to its fiducial unboosted value. We return here to this pr oblem, including 



the p ossibility of a boost factor for the cross section. Our starting point will be IBergstrom et al. 



(|2009l) . where the authors find regions in the (m x , B) parameter space (mass of WIMP vs. boost 
factor) that fits PAMELA/FERMI data based on three different classes of DM models. For sim- 
plicity we will consider only two of those models in our study: from the leptophilic class, the fi 
channel case, where 100% di rect annihilation to U + U~ is assumed and the AHA model, a subclass 



of the Arkani-Hamed model (Arkani-Hamed et al 



2009) which postulates that the Somerfeld en- 



hancement is due to the exchange of a new type of light (sub- GeV) particle. In the AHA case the 
new force carrier is a scalar, (p, which subsequently decays 100% to Somerfeld enhance- 

ment is now generated naturally via ladder diagrams for the XX 00 process, producing 4 muons 
per annihilation versus 2 muons in the direct annihilation case. 
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1.2 Effect of Concentration Parameter 



A second focus of this paper will be to study the dependence of DS properties on the concentration 
parameter of the initial density profile of the halo within which the first stars form. To remind the 
reader, the concentration parameter (c) characterizes how centrally condensed the initial density 
profile is: a larger concentration parameter means that more of the mass is concentrated at the 
center rather than at the outside of the dark matter halo (the precise definition is given below in 
Sec J2.3.2| ). Previous studies of DS considered c = 10 and c = 2. Here we systematically examine 
a sensible range of concentration parameters. The most recent results of numerical simulations of 
structure f ormation seem to favor a "floor" of c = 3 .5 on the concentration parameter of early 



structures ( Zhao et al. 



2003, 



2009 



Tinker et al. 



2010 ) (Some halos can have a smaller concentra- 



tion parameter if star formation begins before the halo is fully formed). At any rate the value will 
differ from halo to halo, even at the same redshift. Thus to explore how the properties of the DS are 
affected by a change in the concentration parameter, we will run the same simulation with three 
different values for it: c = (2, 3.5, 5), which characterizes the overall range for the concentration 
parameter DM halos hosting the first stars. 



1.3 Canonical Values 

We will assume a redshift of formation z = 20 for the first stars, a dark matter halo of 1O 6 M , and 
concentration parameters c = (2, 3.5, 5). 

We will take the annihilation cross section to be 



( av ) = B x 3 x 10" 26 cm7s 
where the boost factor B varies between 100 and 5000 (depend ing on the partic 



(3) 

e phys ics model). 



(2009) described 



The corresponding WIMP mass is taken from the models in iBergstrom et al. 
above and ranges from 100 GeV-4 TeV. 

A key question for DS is the final mass, as the DS accretes more and more material. As long 
as there is a reservoir of DM to heat the DS, the star continues to grow. In the original work 



of 



Spolyar et al.1 (|2009|) . the assumption was made that the initial DM inside the DS annihilates 



awa y in ~ 500 , 000 y ears for a spherical DM halo; here the DS grow to ~ 1000M©. In later work 



of 



Freese et al 



(|2010|) . this assumption was questioned due to the fact that DM haloes are instead 
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triaxial, so that a variety of DM orbits can keep the central DM density higher for longer periods 
of time and the DS can grow supermassive > 10 5 M Q . In reality dark stars will form in a variety 
of dark matter environments and will grow to a variety of masses. For the purpose of illustrating 
how DS vary due to differences in the halo concentration parameter and also due to enhanced 
annihilation rates, we will restrict ourselves to the first option for adiabatic contraction, in which 
the DM originally in the star (due to adiabatic contraction) is the only DM available to the DS; i.e. 
the DS can grow to ~ 1000M Q . 

In addition to this simple adiabatic contraction, we will also consider the effect of captured 
DM on the first stars (Iocco 2008; Freese, Spolyar, and Aguirre 2008). 

In this case, DM passing through the first stars can scatter off the baryons multiple times, lose 
energy, and become bound to the star (direct detection experiments are based on the same physics: 
scattering of DM particles off of nuclei). Subsequently, the star builds up a reservoir made up of 
captured DM, which can power stars. In the 'minimal capture case,' DM heating from captured 
DM and fusion powers the star in equal measure once it reaches the main sequence. 

In section [2] we describe the elements necessary to study the stellar structure of the DS. In 
Section [3] we present results for the influence of varying the concentration parameter and the boost 
factor on the formation and evolution of Dark Stars and their properties. We summarize in Section 

SI 



2 EQUILIBRIUM STRUCTURE OF THE DARK STAR 
2.1 Initial Conditions and Accretion Rates 

The standard picture of Pop. Ill star formation starts with a protostellar 



and cooling via hydrogen cooli ng into a protostar (lOmukai & Nishi 



halo. However, as was found in 



gas c 



oud that is collapsing 



1998) at the centre of the 



Spolyar et al.1 (|2008), DM heating could alter the evolution of the 



first stars significantly. As soon as the protostellar gas reaches a critical core density, DM heating 
dominates over all possible cooling mechanisms. The cloud condenses a bit more but then stop 
collapsing, and becomes a dark star in equilibrium. At this point dark matter annihilation can 
power the dark star. 

As the initial conditions for our simulations we take a DS in which the baryons are fully 
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8 



ionized. This luminous object powered by DM annihilations has a mass of 3M , a radius 1 — 10 
AU, and a central baryon number density ~ 10 17 cm~ 3 . We look for equilibrium solutions as 
described below. From this starting point we follow the ev olution of the DS as it accretes baryonic 



Tan & McKeel (|2004), which decreases 



mass from its surroundings. We use the accretion rate of 
from 1.5 x 10" 2 M o /yr at 3M to 1.5 x lO" 3 M /yr at 1OOOM . At each stage in the accretion 
process we again find solutions in hydrostatic and thermal equilibrium. Eventually the accretion is 
cutoff by feedback effects; as the dark matter runs out due to annihilation, the DS heats up to the 
point where it emits ionizing photons that shut down further accretion. Feedback turns on once the 
surface temperature T e ff reaches 50,000K and is accounted for by introducing a linear reduction 
factor which shuts off accretion completely once the stars surface temperature reaches 10 5 K. 



2.2 Basic Equations 

We u se the numerical c ode previously developed in a paper involving two of the current au- 



thors (ISpolyar et al.H2009Q . Here we review the ingredients in the code before demonstrating the 
modifications relevant to this paper. 

One key requirement is the hydrostatic equilibrium of the star. This is imposed at each time- 
step during the accretion process: 



-P^)—f. (4) 



dP _ f^^M' 

dr r 2 

Above P denotes the pressure, p(r) is the total density and M r is the mass enclosed in a spherical 
shell of radius r. We assume that the star can be described as a polytrope with 



P = Kp 1+1 / n , (5) 



where the "constant" K is determined once we know the total mass and radius (Chandrasekhar 



1939b . 

The energy transport is initially convective with polytropic index n = 3/2 but as the star 
approaches the Zero Age Main Sequence (ZAMS) it becomes radiative with n = 3. The code 
interpolates between n = 3/2 and n = 3 to account for the shift in energy transport as the star 
grows in mass. We can determine the temperature at each point in the radial grid via the equation 
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of state of a gas-radiation mixture, 

= wwtm + i aT(r) ^ p p (6) 

m u n 6 

Here ks is the Boltzmann constant, m u is the atomic mass unit and n = (2X + 3/4F) -1 is the 
mean atomic weight. In all resulting models T ^> 10 4 K except near the surface, so we use the 
mean atomic weight for fully ionized H and He. We will assume a H mass fraction of X = 0.76 
and a He mass fraction Y = 0.24, and that they will remain constant throughout the simulation. 
At each point in the radial grid, T(r) and p{r) a re used to determine the Rosseland mean 



opacity k from a zero metallicity t able from OPAL dig 



low temperatures by opacities from 



Lenzuni et al. 



esias & Rogersl ll996) supplemented at 



(Il99lh for T < 6000K. The location of the 



photosphere is determined by the hydrostatic condition: 

kP = \g (V) 

where g is the acceleration due to gravity at that particular location. This corresponds to a point 
with inward integrated optical depth r ~ 2/3; here the local temperature is set to T e ff and the 
stellar radiated luminosity is therefore: 

L* = ATrRlasT^f, (8) 
with i?* being the photospheric radius. The thermal equilibrium condition is 

L* = Lt t (9) 

where L to t is the total luminosity output from all energy sources as described below in Sec J2.4[ 
Starting with a mass M and an estimate for the outer radius R*, the code integrates Eqns. © 

and © outward from the center. 

The total luminosity output L to t is compared to the stellar radiated luminosity, as in Eq.® and 

the radius is adjusted until the condition of thermal equilibrium is met (a convergence of 1 in 10 4 

is reached). 

2.3 Dark Matter Densities 

2.3.1 Initial Profile and Concentration Parameter 

The first stars form inside ~ 10 6 M Q haloes. Simulations imply that DM halos have a naturally 



cuspy profile, but there is still some uncertainty about the exact inner slope of a DM halo: 
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(2007); 



Springel et al. 



(2008); 



Klypin et al.1 (120101) . Luckily, a previous paper (|Freese et al 



2009) 



showed that a dark star results regardless of the details of the initial density profile, even for the 
extreme case of a cored Burkert profile (such a Burkert pro file is completely u nrealistic). In this 



paper, we use a Navarro, Frenk, & White (NFW) profile (INavarro et al 



1996) for concreteness. 



We assume that initially both the baryons (15% of the mass) and the DM (85% of the mass) can 
be described with the same NFW profile 



p(r) 



(10) 



r/r s (l + r/r s ) 2 ' 

where p is the "central" density and r s is the scale radius. Clearly at any point of the profile, 
baryons will only make up 15% of the mass. The density scale, p Q can be re-expressed in terms of 
the critical density of the universe at a given redshift, p c (z) via 



P0 = Pc{z) 



200 



(ID 



3 /n(l + c)-c/(c+l)' 
where c = r vir /r s is the concentration parameter and r vir is the virial radius of the halo. We will 

assume a flat ACDM universe with current matter density Vt m = 0.24 and dark energy density 

tt A = 0.76. 

One of the main points of this paper is to study the dependence of DS properties on the value 
of the concentration parameter c. Hence we will consider a variety of values, ranging from c=2 to 
c=5. 



2.3.2 Adiabatic Contraction 



As the baryons start to collapse into a protostellar cloud at the center of the DM halo, the DM 
responds to the changing gravitational potential well and becomes compressed. As described in 
our previous work, we will use adiabatic contraction (AC) to calculate the effect of baryons on the 
DM profile. With an initial DM and gas profile and a final gas profile, w e can use the adiabatic 



invariants to solve for the final DM profi 



Blumenthal et al. 



1986 



Rvden & Gunn 



e. We use the Blumenthal method (Barnes & White 



19841 : 



1987|) to calculate the adiabatic compression of the halo. 
In this case the simplifying assumption of circular orbits is made. Angular momentum is the only 
non-zero invariant. Its conservation implies that M/(r/)r/ = Mj(rj)?Y In the case of circular 
orbits M is the mass interior to the radius r of an orbit and the indices / and i refer to final and 
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initial orbits respectively. As mass grows inside of the orbit, its radius must shrink and the DM 
profile steepens. 



Freeseetal 



(2009), where a 



The validity of this method in this conte xt has b een checked in 
more precise algorithm, developed by Young ( Youngtfl98oL has been used and a difference within 
at most a factor of 2 has been found. Whereas the Blumenthal method assumes circular orbits, 
Young's method only assumes spherical symmetry of the system. Therefore only one of the three 
conserved actions is identically zero in this case. Namely the plane of each orbit does not change. 
The other two action s, angular moment um and the radial action respectively are non-zero and 



conserved. In view of 



Freese et al. 



(|2009D we are confident that the simple Blumenthal method is 
sufficiently accurate for our purpose. 

In this paper, we assume that all of the DM moves on circular orbits. The DM will become 
exhausted once all of the DM on orbits interior to the DS have been depleted. The timescale for 
this to happen is on the order a mil lion years. This is probably an unduly cautious assumption. 
In a recent paper ( Freese et al. 201ol) we studied the case of triaxial haloes with large numbers of 
centrophilic orbits, which provided a much larger DM reservoir than we are considering, leading 
to supermassive DS. In reality dark stars will form in a variety of dark matter environments and 
will grow to a variety of masses. In this paper we could have grown the stars to supermassive sizes, 
but do not consider it necessary for showing the dependence of DS properties on the two effects 
we are studying. 



2.3.3 Dark Matter Capture 



Until now we have only discussed the DM inside the DS due to adiabatic contraction. However, 
dark star's DM reservoir can be refueled by DM capture. This refueling requires an additional 
piece of particle physics: scattering of DM off the atomic nuclei inside the star. Some of the 
WIMPs from far out in the halo have orbits passing through the star. This DM can be captured 
and sink to the stars center, where they can annihilate efficiently. The capture process is irrelevant 
during the early evolutionary stage of the DS, since the baryon density is not high enough to 
efficiently capture DM. However, at the later stages as the DS contracts towards the ZAMS, the 
baryon densities become high e nough for su bstan tial capture to be po ssible. This mechanism was 



first noticed simultaneously by (|locco 



2008) and ( Freese et al 



2008b). 
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The capture rate is sensitive to two uncertain quantities: the scattering cross section of WIMP 
interactions with the nuclei and the background DM density. In terms of the relevant particle 
physics, we will consider only spin-dependent (SD) scattering cross sections with 

a c = l(T 39 cm 2 . (12) 

The details of our capture study have previously been presented in Freese et al, 



(2008b) and 



will not be repeated here. We do wis h to emphasize that w e will assume the same case of minimal 
capture that we previously studied in lSpolyar et al.l (I2009|) . in which case the heating from fusion 
and from DM heating are taken to be comparable. The more extreme and interesting case of 
dominant capture, which could last as long as the DS continues to exist in a high DM density 



environment, was studied elsewhere (IFreese et al, 
DS. 



2010) and can lead to supermassive > 1O 5 M 



2.4 Energy Sources 

There are four possible contributions to the DS luminosity: 

grav * -^nuc ~ -^cap (13) 

from DM annihilation, gravitational contraction, nuclear fusion, and captured DM respectively. 
2.4.1 DM Annihilation 

The heating due to DM annihilation is given in Eqns. (QQ) and © and dominates from the time of 
DS formation until the adiabatically contracted DM runs out. In order to compute the luminosity 
generated by DM heating, one needs to know what fraction of the total energy generated by WIMP 
annihilations is deposited in the star. This quantity, which we n ame /q, will be different for various 



models for DM. In previous work (e.g. in lSpolyar et al.l (I2009|) ) a fiducial value for fg of 2/3 was 



used as appropriate for many typical WIMPS, under the following assumptions. First, the final 
products of DM annihilation, after all unstable particles have decayed to the lightest states, are 
taken to be three types: i) stable charged particles (i.e. e*), ii) photons and iii) neutrinos. Second, 
the energy distribution was assumed to be roughly equal among the three final products listed 
above. The electrons and photons were taken to thermalize with the star and deposit their energy, 
whereas the neutrinos escape; hence fg ~ 2/3. 
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On the other hand, given a specific model for DM, one could compute the precise value for fg 
using the energy distribution of all final annihilation products for the model under consideration. 
While this procedure could be done for the specific models in this paper, still we will use the 
standard /q = 2/3 in order to make simple comparisons with our previous work[]. Since the 
aim of this paper is to understand the ef fect of the boost f actor on the DS properties, we will fix 



the value of /q to 2/3, the same as in ISpolyar et al. 



(2009). The differences between the boosted 



and unboosted cases will now be due only to the boost factor itself and the different masses of 
the WIMP in the two cases (rather than to the detailed values of /q). Since /q appears always 
multiplied by (av ), any ambiguity we have introduced by fixing the energy deposition factor to 
2/3 can be traded for an ambiguity in the annihilation cross section. 



2.4.2 Gravitational energy 



Once the DS runs out of DM, it begins to contract; gravity turns on and powers the star. Although 
relatively short, this Kelvin-Helmholtz (KH) contraction phase has important consequences: it 
drives up the baryon density and increases the temperature, therefore leading to an environment 
where nuclear fusion can take place. For a polytr ope of index n = 3 the gravitational contribution 

(|2009r) using the virial theorem, 

3 d ( GM 2 ' 



to the luminosity was found in 



Freese et al 



J grav 



Adt 



R 



--E d 

2dt radl 



(14) 



where E rad = J dVaT 4 is the radiation energy. 



2.4.3 Nuclear Reactions 

Subsequently nuclear reactions become important. We include the following nuclear reactions, 
which are typical for a zero metallicity star during the pre main sequence evolution: i- burning of 
primordial deuterium (assumed to have a mass fraction of 4 x 1CT 5 ) which turns on rapidly once 
the stars central temperature reaches ~ 10 6 K, ii- the equilibrium proton-proton cycle for hydrogen 
burning, and Hi- the alpha-alpha reaction for helium burning. Since we track the evolution of the 

3 For the leptophilic models we are considering in this paper, DM annihilates either directly or via exchange of a light scalar to so that 

one might more reasonably expect smaller values since the muons decay to electrons and two neutrinos. The actual value will be closer to ~ 1/3 
but will depend on the energy distribution among the various final products; still it will be a number of 0(1). 
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DS only until it reaches ZAMS we do not need to consider the changes in stars atomic abundances. 
In order to calculate the nuclear luminosity, defined as 



we use the methods described in 



dMe r 



(15) 



Claytonl (|1968|) to obtain the energy gene ration rate, e nur .- For the 



proton-proton reaction we u se the astrophysical cros s sectio n factors from 
He burning parameters from lKippenhahn & Weigertl (|l990|) . 



Bahcall 



(1989) and the 



2.4.4 Luminosity due to captured DM 

As we have seen in Sec J2.3.3l during the later stages of the pre main sequence evolution captured 
DM can become an important energy source. The luminosity due to DM capture is 

Leap ^Tfl^T ca p, (16) 

where 

T cap = f Q J dVp 2 cap (av)/m x (17) 

and the factor of 2 in Eq.(fT6l) reflects the fact that the energy per annihilation is twice the WIMP 
mass. In all simulations we will consider the case of "minimal capture", which corresponds to 
equal contribution to the luminosity from capture and nuclear fusion when the star reaches the 
main sequence. 



3 RESULTS 

On the whole, for all the values of boost factor and concentration parameter we have considered 
in this paper, the results are roughly the same: the final DS is roughly ~ 1OOOM , ~ 1O 7 L , and 
lives ~ 10 6 yrs. Thus if the e + excess seen in PAMELA is due to WIMP annihilation, the required 
leptophilic boosted cross section is certainly compatible with the DS picture. However there are 
interesting differences between models which we will discuss. 

Other than in the subsection immediately following this one, we will consider four WIMP 
models. As motivated below, we will focus on one boosted model denoted by AH4 with the fol- 
lowing set of parameters: B = 1500, m x = 2.35 TeV and c = 3.5. As our unboosted models, we 
will take 100 GeV WIMPs with the canonical cross section of 3 x 10~ 26 cm 3 /s , and we will con- 
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sider three values of the concentr a tion p arameter, c = (2, 3.5, 5). For comparison, the "canonical 



case" considered in 



Spolyar et al. 



(2009) was the unboosted 100 GeV case with c=2. The relative 
boost factor between the AH4 model and the unboosted models is best described as follows: since 
Eq.(l) tells us that DM heating scales as (av) /M x , one can see that the AH4 is exactly equivalent 
to a 100 GeV WIMP mass with Boost Factor 150/2.35 ~ 64. In other words the relative boost 
factor between the AH4 model and the unboosted cases is actually 64. 

Also, we assume throughout that the boost factor will have the same value for (av) that is 
required to explain the PAMELA/ Fermi data. Particle physics models written down can certainly 
allow for larger boost factors in the halos in which the first stars form (the converse does not hold; 
DM has a lower velocity in a DS halo versus the Milky Way halo, which would act to increase 
the Sommerfeld enh ancement). For clarity, we will fix the boost factor to the values given by 



Bergstrom et al. 



(2009). 



3.1 'Final' Luminosity for Dark Stars with boosted DM 



We will investigate the 'final' luminosities and masses of Dark Stars as they enter the main se- 
quence. We study a variety of WIM P parameter ranges c apable of explaining the PAMELA data; 



in particular we follow the work of iBergstrom et al 



(2009) and take various boosted DM annihila- 



tion rates and WIMP masses from fig. 1 or their paper, which gives 2a contours in the enhancement 
factor - mass plane needed to fit PAMELA and Fermi data. In this section we consider two of the 
particle physics models they consider: the fi channel and AH4 type of models. 

In Fig. Q] we plot the 'final' luminosity at the end of the DS lifetime, when the DS has ac- 
creted to its maximum mass: the dark matter has run out and the star is about to enter the zero 
age main sequence where it will be powered by fusion. We plot the luminosity as a function of 
the boost factor. The points represent individual simulation outputs. The solid (dashed) lines con- 
nect simul ation outputs generated with the DM mass and Boost factors from the 2a contours in 



fig.l from 



Bergstrom et al. 



Berg strom et al.l (|2009h corresponding to the Fermi (PAMELA) confidence regions in 



(|2009|) . We note that even if the boost factors range over more than one order 
of magnitude, from 100 to 5000 ( the corresponding WIMP masses take values in the 1-4 TeV 
range) the final luminosities in all cases are relatively similar, ranging from ~ 7 x 10 6 L Q to 
~ 9 x 1O 6 L (The variation in luminosity for a fixed boost factor corresponds to the range of 

© 0000 RAS, MNRAS 000, 000-000 



16 



1 x10 7 



9x10 6 



w 8x10 

o 

c 

E 



7x10 e 



6x10 e 



fji channel mode 
i 



AH4 model 




1 x10 7 



100 



- 9xio e 



- 8x10° 



- 7x10 e 



J6x10 e 



1000 
Boost 



100 



1000 



Boost 



Figur e 1. "Final" luminosities (as the Dark Star enters the ZAMS) as a function of Boost factor for the /^-channel and AH4 models examined 
in(] 
2o 

model and is consistent with both Fermi and PAMELA. 



rigure i. rinai luminosities ^as ine uarK oiar eniers ine z^rtivi^j as a lunction oi Boost lactor lor tne /i-cnannei anu /\n^ moueis examineu 
in Berg strom et all 12009]). The contour regions are generated using the range of WIMP masses for a fixed boost factor taken from the corresponding 
2a contours in fig. 1 in lBergstrom et alj 12009). The central star-shaped point in the right panel will be taken to be our designated boosted AH4 



allowed m x in iBergstrom et al.1 (120091) . For instance for a boost factor of 1000, the upper bound 
of the Fermi contour in left panel of Fig. Q] corresponds to a m x = 1.4 TeV; the lower bound 
to m x = 1.9 TeV ). As the boost factor increases, we note that the final luminosity gets slightly 
smaller. This is due to the fact that the ratio (av) /Mdm typically also increases in most cases. This 
leads to a faster depletion of the adiabatically contracted DM from the star, shortening the time un- 
til it reaches the main sequence, therefore reducing the amount of mass accreted and consequently 
the luminosity at that point. 
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3.2 Structure and Evolution of DS 

In this section, we will analyze the two effects of i) boosting the annihilation cross section and ii) 
a variety of concentration parameters on the structure of a Dark Star. As mentioned above, for the 
"boosted" case we choose one sample point from the AH4 type models, which corresponds to the 
large star shaped point at the center of the right panel of Fig. [Q Henceforth we will denote by AH4 
this point which has B = 1500, m x = 2.35 TeV and c = 3.5. 



3.2.1 Luminosity Evolution 

In Fig. [2] we compare the luminosity evolution for the four cases we have just described. In all 
cases, DM annihilation heating provides the dominant contribution to the DS luminosity until the 
DM runs out. At this point Kelvin-Helmholtz contraction sets in and briefly provides the dominant 
heat source, until the star becomes hot and dense enough for fusion to begin. In the cases where 
we include capture, DM annihilation may again become important at the later stages leading to a 
new DS phase. 

In both upper panels of Fig. [2] we take c = 3.5: the left panel is the boosted AH4 model while 
the right panel is unboosted (the relative ratio of boost factors is effectively 64 as mentioned 
before). Due to the ambiguity in the value of the concentration parameter we will also study 
its implications on the evolution of the Dark Star by running the unboosted 100 GeV case for 
c = (2, 3.5, 5) respectively; the lower two panels are the un boosted case with c = 2 and 5. The 



lower left panel is the same as the canonical case studied in lSpolyar et al. 



(2009). 



Effects of Boost: The DM heating is the most powerful in the boosted AH4 case, since Eq.(l) 
indicates that heating scales with cross section. Thus at any given time during the DS phase, this 
model has the brightest luminosity, as can be seen in Fig. [21 Consequently, the DM is burned up 
the more quickly in the boosted AH4 case, leading to the shortest DS lifetime. 

In order to balance the higher DM heating, a larger radius is required, which leads to a lower 
central temperature and density. See Fig. [9] As the adiabatically contracted DM runs out the evo- 
lution is similar in all cases, the final luminosity being on the order of 10 7 L Q . 

Effects of concentration parameter: When we compare the three different concentration pa- 
rameter cases, we notice that the nuclear fusion sets in earlier for smaller values of c. The higher 

© 0000 RAS, MNRAS 000, 000-000 



18 



10 



10 



10 E 



10 e 



10* 



10 2 



10 c 



AH4 



totd 

grovity 

nuclear 

capture 

annihll. 



c = 3.5 




T 7 



totol # 

gravity 

nuclear 

capture 

annihll. 




10 



10 B 



10 e 



10 4 



10 2 



10° 



10 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 
c = 2.0 c = 5.0 



10 



10 



10 E 



total 

gravity 

nuclear 

capture 

annihll. 




total 

gravity 

nuclear 

capture 

annihll. 




i 
i 
i 

■ i L 



10 



10 8 



10 6 



10 4 



10 2 



10° 



10 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 
t [Myr] t [Myr] 

Figure 2. Luminosity evolution as a function of time for the four cases under consideration. The upper left panel displays the only boosted case, 
the AH4 model, with B = 1500, m x = 2.35 TeV and c = 3.5. The other panels correspond to the unboosted case of a 100 GeV WIMP with 
(<tv) = 3 X 10 26 cm 3 /s and concentration parameters as labeled. 
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the concentration parameter the more adiabatically contracted DM available, therefore it takes 
longer to start the transition to the main sequence. This will lead to slightly higher final masses, 
which in turn translate to higher final luminosities. Nevertheless, as pointed out before they are all 
very close to 1O 7 L . In all three cases where (av) = 3 x 10" 26 cm 3 / s we notice a flash in the lumi- 
nosity when both the gravitational energy and the energy due to annihilations of the adiabatically 
contracted DM are relevant. This happens at a time somewhere between 0.31 Myr (for c = 2 ) and 
0.42 Myr (for c = 5). The same flash can be seen for the AH4 panel, but now the true maximum 
of the luminosity occurs after only 0.1 Myr. The maximum luminosity is ~ 4 x 1O 9 L . 

The main differences between the different WIMP models as regards the luminosity evolution 
in Figf2]are in the time the "pure" dark star phase lasts. The higher the boost factor, the shorter 
this phase. Conversely a larger concentration parameter c prolongs the DS phase, since more DM 
is available. 

3.2.2 Pre main sequence evolution, HR diagram 

In Fig. [3] we plot Hertzsprung-Russell (H-R) diagrams for the four cases. One can see two distinct 
phases. First, the DS goes up the Hayashi track with a very steep increase of the luminosity yet 
relatively cool surface temperature, T e ff ^ 10 4 K. At the end of the Hayashi track the star enters 
the Henyey track. This path corresponds to an almost constant luminosity while the temperature 
increases fast, mostly due to the KH contraction phase. As a rule of thumb once a star is on the 
Heyney track its core should be fully radiative. The graphs end at a temperature of ~ 10 5 K when 
the star reaches the main sequence. 

Effects of Boost: From the left panel of Fig. [3] one can see that the boosted AH4 case has the 
highest luminosity, due to the extremely efficient DM heating Qdm ~ (&v) /m x forming a lumi- 
nosity peak. However, as the AH4 case burns up its DM, the its luminosity falls. The boosted and 
unboosted cases eventually cross over at a temperature of ~ 10 4 K, and henceforth the unboosted 
case has a higher luminosity. Consequently, the boosted AH4 case has the lowest luminosity as the 
star moves onto the main sequence, as discussed above. 

Effects of concentration parameter: In the right panel of Fig.[3]the trend is uniform: an increase 
in the concentration parameter leads to an increase in the luminosity. The difference is relatively 
small in the early stages of the evolution, at low temperatures. This is due to the fact that the 
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Figure 3. Hertzsprung-Russell diagram for DS. The left panel displays the unboosted 100 GeV case (dashed line) and the AH4 case, both for 
c=3.5. The right panel displays the unboosted 100 GeV case for a variety of concentration parameters as labeled. Also labeled are a series of points 
corresponding to the evolution of the DS towars the ZAMS at different times. 

adiabatically contracted DM density profile is not very sensitive to the concentration parameter, 
therefore about the same amount of DM heating will be generated in each case. However, for a 
lower value of the concentration parameter, the adiabatically contracted DM runs out faster, as 
there is less DM available. This leads to a shorter "pure" DS phase, as can also be seen from Fig. 
[2l and consequently to slightly lower final mass and luminosities. 



3.2.3 Radius and Effective Temperature ofDS 

Effects of Boost on Radii: In the "pure" DS phase the star will be puffier for a higher (av)/m x 
ratio, as can be seen from the left panel of FigH This is expected due to the much higher DM 
heating in that case: a larger radius is required to balance the DM energy production and the 
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Figure 4. DS radius as a function of time. The left panel displays the unboosted 100 GeV case (dashed line) and the AH4 case (solid line), both 
for c = 3.5 The right panel corresponds to various concentration parameters c, for the unboosted 100 GeV WIMP. 



radiated luminosity, which scales as R 2 . For instance, in the boosted AH4 case the maximum 
radius is at about 2 x 10 14 cm, whereas for the 100 GeV non-boosted WIMP, the DS will have a 
maximum radius of 10 14 cm. However, as mentioned before, the KH contraction will set in earlier 
in the boosted case. This phase corresponds to the sharp decrease in radius in FigJH The final radii, 
as the DS enters the ZAMS are similar in both cases, at around 6 x 10 11 cm. 

Effects of Concentration Parameter on Radii: From the right panel of FigJH we notice an uni- 
form increase in the radius with the concentration parameter. Again, more efficient DM heating 
leads to a puffier DS. The maximum radii range from 8 x 10 13 cm to 10 14 cm for c = (2 — 5). After 
the KH contraction phase the DS settles to a radius very close to 6 x 10 11 cm as it enters the ZAMS. 
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Effective Temperatures: In Fig. \5\ we plot the effective surface temperature as a function of 
time. At first, during the DS phase, T e ff is relatively constant below 10 4 K. Once the DM starts 
to run out, the star contracts and heats up, leading to a sharp increase in T e ff due to the onset of 
the KH contraction phase. Once nuclear fusion becomes the dominant energy supply and the star 
ceases to contract, the surface temperature reaches a plateau. The final value of T e ff ~ 10 5 K is 
always the same for all cases, regardless of value of boost factor or concentration parameter. The 
AH4 case (left panel in FigJ5]) starts as a cooler star, again typical of more efficient DM heating 
in that case due to the boost factor. In addition, in the AH4 case the DM runs out more quickly, 
leading to an earlier increase in temperature. Regarding the various concentration parameters: 
during the DS phase, the surface temperature is roughly the same in all cases, ~ 7 x 10 3 K. The 
DS phase lasts the longest for the highest value of c as this case has the most DM to begin with; 
thus the temperature starts to rise later for ever larger concentration parameters c. 

3.2.4 Energy Transport near the Core 

The DS starts with a fully convective structure, modeled by a fluid with a poly tropic index n = 1.5. 
Then a radiative core starts to develop, that grows outwards, until most of the star is described by a 
polytrope of index n = 3. In Fig.[6]we plot the radiative gradient in the innermost zone at the center 
of the DS. The dashed horizontal line illustrates the critical value for convection. Models above the 
line have a convective core, while models below the line have radiative cores. Models with more 
efficient DM heating - i.e. the models with higher values of c or the AH4 model — transition to 
radiative energy transport later; stars with more efficient DM heating require a larger radius. At a 
fixed luminosity with more efficient DM heating, the star must have a larger radius and a lower 
densities to keep DM heating and the stellar luminosity in balance. With a lower density, the star 
has a lower central temperatures. At lower temperatures, the number of bound states increases, 
which increases the number of bound-free transitions. Also, the number of free-free transitions 
increase. These two effects increase the opacity which produces a larger radiative gradient and 
delays the transition from convective to radiative transport. 

Once the original DM in the star runs out and nuclear fusion begins, a convective core develops 
in all cases. The central gradient is high due to the fact that nuclear fusion takes place primarily 
in the core of the star. Similarly, once the star repopulates its DM due to capture, this new DM 
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Figure 5. Effective temperature as a function of time. The left panel corresponds to the unboosted 100 GeV case (dashed line) and the boosted 
AH4 WIMP parameters defined previously (the solid line), both for c = 3.5. The right panel corresponds to various concentration parameters c, 
for the unboosted 100 GeV WIMP. 

population is thermalized with the star and its density is also sharply peaked at the center of the 
star. Thus both nuclear fusion and captured DM lead to a large radiative gradient in the core and 
therefore favor convection. 

Again, the transition happens later for higher values of c. The time when the convective core 
develops corresponds almost precisely to the time when the captured DM heating becomes signif- 
icant, as can be seen by comparing with Fig. [21 0.31 Myr, 0.37 Myr and 0.42 Myr for c = 2, 3.5 
and 5 respectively. 



3.2.5 Baryonic Central Density 

The baryon central density is plotted in Fig.[7]for the four cases we considered. 
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Figure 6. Radiative gradient at the center of the DS as a function of time. Models with the gradient above the horizontal line are unstable to 
convection. The lines labeled by the value of the concentration parameter c correspond to a DS powered by 100 GeV WIMPs with canonical 
unboosted cross section. The model AH4 is as defined previously in Fig. 1. 



The higher cross section of the AH4 case at first leads to a puffier star (larger radius to keep 
the radiated luminosity at a level to balance the higher DM heating); thus it is not surprising that 
the AH4 case initially has a relatively lower central density. However, the initial DM in the DS 
runs out earlier in the AH4 case due to more efficient burning; hence in the left panel the two lines 
cross due to the earlier onset of the KH contraction (marked in the plot by the sharp increase of 
the densities) in the AH4 case. 

The central density p&(0) scales inversely with the concentration parameter, as can be seen 
from the right panel of the same plot. Again, more DM heating (higher c) will lead to larger radii, 
therefore smaller densities. However, as opposed to the situation depicted in the left panel, In the 
right panel, the curves do not cross since the models with a larger concentration parameter have 
more DM, which delays the onset Helmholtz contraction. We will come back to this in Sec J3.2.7l 
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Figure 7. Central baryon density as a function of time. The left panel displays the unboosted 100 GeV WIMP case (dashed line) and the boosted 
AH4 case (the solid line), both for c=3.5. The right panel corresponds to various concentration parameters c for the 100 GeV WIMP case with 
unboosted canonical cross section. 

In all cases considered here once the star goes onto the main sequence, the central density is close 
to 100g/cm 3 . 

3.2.6 Mass as DS enters the ZAMS 

In Fig. [81 we have plotted the DS mass as a function of time. In all cases, the final mass when 
the DS enters the main sequence is ~ 700 — 1OOOM ; however, there are slight differences for 
different models. The models with more effective DM heating - i.e. the AH4 model compared to 
the 100 GeV unboosted case - burn up their original adiabatically contracted DM the most quickly 
and enter the KH contraction phase the soonest. This in turn leads to a smaller final mass, as feed- 
back effects will shut off accretion sooner. When comparing the cases with different concentration 
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Figure 8. DS mass as a function of time. The left panel displays the unboosted 100 GeV case (dashed line) and the AH4 model, both for c=3.5. 
The right panel displays the unboosted 100 GeV case for a variety of concentration parameters c. 



parameter we notice that the DS final mass is an increasing function of c. As previously explained 
an increase in the concentration parameter leads to a longer DS phase, and hence to more mass 
accreted. For the case of the unboosted 100 GeV WIMP with c = 3.5 the final mass is around 
900M Q , whereas for the AHA case it is close to 7OOM . For c ranging from 2 — 5 the DS will have 
a mass in the 800M Q — 1OOOM as it reaches the main sequence. 



3.2. 7 Density Profiles for DM and Baryons inside the DS; amount of adiabatically contracted 
DM 

In Fig. [9] and FigJTOl we have plotted the density profiles of the adiabatically contracted DM and 
baryons respectively. The AH4 model for the same stellar mass has a lower DM density than 
the canonical unboosted 100 GeV case by roughly an order of magnitude, and also has a more 
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extended profile. For instance, in the case of a DS of 3OOM the values are 5 x 10 _11 g/cm 3 (AH4) 
and ~ 1 x lCT 9 g/cm 3 (canonical) respectively. This effect can be attributed to the fact that at a 
higher annihilation cross sections or a lower particle massesj, a larger radius and a lower density 
are needed to balance the DM heating and the stellar luminosity (which scales as R 2 ). 

During the early stages of the DS evolution the dependence of the adiabatically contracted DM 
density on the concentration parameter is very small at least for the range we have considered here. 
As it can be seen from FigJTl prior to the onset of the KH contraction phase, the central baryon 
densities for models with different values of the concentration parameter have similar baryon den- 
sity and DM density as well. Models with a larger concentration parameter have slightly more dark 
matter, have slightly lower central DM densities, and are also more extended. Before the contrac- 
tion phase, the central DM density of the c = 5 case's density is 10% lower than the c = 2 case. 
The radius is 20% larger. Models with different concentration parameters only begin to dramati- 
cally diverge once the star begins to contract and enters its Kelvin-Helmholz contraction phase. At 
this point, the star begins to shrink, which cause the DM densities to increase dramatically. 

DS in halos with a larger concentration parameter have more DM and thus delay the onset of 
the KH phase. For c = 2, the contraction phase begins at t ~ 0.28 Myr; See FigJU At this time, 
the stellar mass has reached ~ 7OOM . For c = 3.5 and 5, the contraction phase begins later. In 
the case of c = 5, the star has a mass of ~ 850M Q . Thus the contraction begins once the star is 
more massive. 

At a fixed stellar mass, the DM densities will differ dramatically between models which are in 
the contraction phase compared to those which are not contacting . For instance, let us consider a 
750M o DS. In the case of c = 2, the star has entered the contraction phase. While the cases with 
a larger concentration parameter (3.5,5) have yet to begin their contraction phase. Hence the stars 
with a larger concentration, have a lower DM density and are more extended, which can be seen 
in Fig.|U 

Finally, in Fig. [TT|we have plotted the amount of adiabatically contracted DM inside the DS as 
a function of time. One can also see that DM densities are many orders of magnitude lower than 
their baryonic counterparts at all times. Although the amount of DM never exceeds 0.4M Q , yet 



4 A similar effect was noticed in Spolvar et al. 2009), where it was found that "the average DM density in the star is an increasing function of 
M x "; n.b. the higher annihilation cross section of the AH4 case can be traded for a lower WIMP mass since heating scales as (av)/M x . 
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Figure 9. Adiabatically contracted DM density profiles. Each line corresponds to a fixed value of the mass of the DS during its evolution. Note that 
certain lines that are mentioned in the legend do not appear plotted in all four panels. This is due to the fact that at that stage the DS has exhausted 
all the DM. 
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Figure 10. Baryonic density profiles at different stellar masses for the four cases under consideration. The solid line corresponds to the mass 
DS as it enters the ZAMS. 
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this is sufficient to power the DS all the way up to ~ 1OOOM (where most of the mass is baryons) 
and 1O 7 L . Indeed DM heating is a very powerful energy source. 
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4 SUMMARY AND CONCLUSIONS 

In this paper we have considered the effect on DS of leptophilic boosted DM annihilation cross 
sections, as would typically be required to explain PAMELA data. Second, we have varied the 
concentration parameter in a host of DS models. We have restricted our study to include the fol- 
lowing two sources of DM: i) the DM originally contained in the star due to adiabatic contraction 
and ii) the minimal capture scenario. We have not included the additional DM du e to extended 
adiabatic contraction or to maximal capture models discussed in lFreese et al.1 (|2010|) . Nonetheless, 
the dependence of DS properties on boost factor and concentration parameter can easily be seen in 
this paper. As our prototypical boosted case, we have focused on the AH4 model with the follow- 
ing parameters: B = 1500, m x = 2.35 TeV and c = 3.5. In the unboosted cases, we have taken 
M x = 100 GeV and three values of the concentration parameter, c = 2, 3.5, and 5. 

We have found that, if the positron excess observed by PAMELA is indeed due to leptophilic 
boosted DM, then there would be an early DS phase of stellar evolution powered by DM heating, 
lasting long enough to bring the star to substantially higher mass and luminosity than predicted 
for regular Pop.III zero metallicity stars. Our basic results are that the final stellar properties, after 
the DS runs out of or its original adiabatically contracted DM fuel, undergoes Kelvin Helmholtz 
contraction and enters the main sequence, are always roughly the same: ~ 1000M o , ~ 10 7 L Q , 
lifetime ~ 10 6 yrs. Similarly, these same DS properties are also the basic result independent of the 
value of the concentration parameter in the range between c=2 and c=5. 

We reiterate that these values are only for the case of simple adiabatic contraction and minimal 
capture; if we were to include the additional DM due to extended AC or maximal capture, then 
these values would be different by many orders of magnitude. However, the basic dependence 
on the parameters of interest in this paper would generalize. In particular, the result that the final 
mass, luminosity, and lifetime are relatively independent of boosted cross section or value of c 
would still hold up. In addition, the slight differences from model to model, discussed in the next 
paragraph, would also hold up. 

We have found that the lifetime, final mass, and final luminosity of the DS, though roughly 
similar in all cases, do show some dependence on boost factor and concentration parameter. We 
have found that the DS lifetime is shorter in the boosted case, since the DM is exhausted sooner. 
On the other hand, the lifetime is longer for higher concentration parameter since there is more 
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DM to begin with. Thus nuclear burning starts earlier in the case of unboosted cross section or 
low concentration parameter. The DS accretes matter continuously while it remains powered by 
DM heating. Hence the largest final stellar mass results for the longest living DS, i.e. the un- 
boosted case with the highest values of c. In all cases, the final mass is ~ 1OOOM . We have 
shown the H-R diagram for all cases, studying both the Hayashi track and the approach to fusion 
power. We have showed that the 'final' luminosity at the end of the DS lifetime is ~ 1O 7 L in 
all cases, with the luminosity decreasing slightly as a function of increasing boost factor or de- 
creasing c, again because of the more rapidly depleted pool of DM. We have also examined the 
luminosity evolution of the DS as a function of time. During the DS phase itself, the luminosity is 
higher/lower for boosted/unboosted cross sections. The reduced luminosity during the DS evolu- 
tion in the unboosted case is a consequence of the reduced energy production, since DM heating 
is proportional to annihilation cross section and the square of the DM density. Then at lower cross 
section (unboosted), a smaller radius is needed to balance the DM energy production and the ra- 
diated luminosity. Consequently the unboosted case has higher central densities (both for baryons 
and DM). Similarly, lower values of c have lower luminosity during the DS phase, smaller radii, 
and higher central density. In all cases the DM density is a minute fraction of the total density, 
with baryons dominating the gravitational potential; we have shown the density profiles of both 
components. Again in all cases, the amount of DM inside the star never amounts to more than 
O.4M , a tiny fraction of a star that grows to ~ 1000M Q ; yet this DM is sufficient to power the 
star. This is "the power of darkness." 
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